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Abstract. The TeV-emitting BL Lac object Mkn 421 was observed with very long baseline interferometry (VLBI) at three 
closely-spaced epochs one-month apart in March-April 1998. The source was also monitored at very-high y-ray energies 
(TeV measurements) during the same period in an attempt to search for correlations between TeV variability and the evolution 
of the radio morphology on parsec scales. While the VLBI maps show no temporal changes in the Mkn 421 VLBI jet, there is 
strong evidence of complex variability in both the total and polarized fluxes of the VLBI core of Mkn 421 and in its spectrum 
over the two-month span of our data. The high-energy measurements indicate that the overall TeV activity of the source was 
rising during this period, with a y-ray flare detected just three days prior to our second VLBI observing run. Although no firm 
correlation can be established, our data suggest that the two phenomena (TeV activity and VLBI core variability) are connected, 
with the VLBI core at 22 GHz being the self-absorbed radio counterpart of synchrotron self-Compton (SSC) emission at high 
energies. Based on the size of the VLBI core, we could derive an upper limit of 0.1 pc (3 x 10 17 cm) for the projected size of 
the SSC zone. This determination is the first model-free estimate of the size of the y-ray emitting region in a blazar. 
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1. Introduction 

Among active galactic nuclei (AGN), BL Lac objects form 
a class of sources characterized by high radio and opti- 
cal variability, dominance of continuum over line emission, 
and strong and variable polarization. Such extreme proper- 
ties have already suggested long ago that substantial rela- 
tivistic beaming most pro b ably occurs in this type of source 
(Bla ndford&Reesl Il978t iBlandford & Konigll Il979i) . This 
idea was confirmed by direct detection of apparent super- 
luminal motion in the radio jets of many BL Lac objects 
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1998; 


Fossati et alJl998: Sambruna et al. 2000: Bai & Leel2001 ). 



line interferometry (VLBI) technique ( Phillips & Mute] 


1982 


Witzel et al.| [l988; 


Gabuzda et al. 


1989; 


Mutel et alJ 


1990 


Gabuzda et al.ll200C 


. Building on these findings, AGN unifi- 



cation schemes have emerged, all basically describing BL Lac 
objects as radio-loud AGN with relativ istic jets pointing within 
a small angle towards the observer dUrrv & Padovaniljl995 ; 
iFalckeet all 19911 lGopal-Krishnall995UKollgaard et all l 996 ; 



In the past decade, more evidence of relativistic beaming 
came from high-energy X-ray and y-ray data, which revealed 
intense fluxes and strong variability up to the TeV level, 
especially in the two nearby BL Lac objects Mkn 421 and 
Mkn 501 dPunch et alJl992HOuinn et al. Il9 96). In the case of 
Mkn 421, very high relativistic Doppler factors (between 10 
and 15) are demanded in order to reproduce the dramatic rapid 
flares that have been observed on timescales as short as 3 min 
dGaidos et alJ 1 19961 ICoppil 1 19971 ICelotti et alJ Il998l Icul 
20041 lAharonian et alJ l2005). Similar or even higher Doppler 
factors (in the range 20-50) are also obtained when fitting 
basic synchrotron self-Compto n models to the b roadband 
spectr u m of Mkn 421 fMaraschi et alJll999t iTakaha shi et alJ 
200 0l iKrawczvnski et alJ 1200 it iKonopelko et alJ l2003t 



Katarzvnski et all2003l) ~ 
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The VLBI observations have sought evidence of superlu- 
minal motion in Mkn 421 ever since the early 1980's. A first 
series of VLBI maps at 5 GHz reported apparent motions of 
about 2c and an angle to the line of sight of 34° for the parsec - 
scale iet (Baat h et alJl98lllBaatlJl984llzh"a ng & Baath 1990). 
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On the other hand, the source was found to be unresolved at 
high frequ ency with a core size of 0.15 milliarcsecond (mas) 
at 22 GHz (|zhang & B aathll99ll) . In the 1990's, the improved 
performances of the VLBI technique p ermitted the detection 
of a weak one-sided iet JPolatidis et alJll995l: Iku et alJll99"l 
Edw ards et al.lll998t iKellermann et alJll998t iGiovannini et alJ 
19991) . This jet shows wiggles starting at about 5 mas from the 



core, as well as strong distorti ons at a distance o f 20 mas from 
the core. Based on these data. IGiovannini et alJ dl999h derived 
a viewing angle smaller than 30° and an apparent jet speed be- 
tween ~0.8c and lc. While faster speeds (~2c) were reported, 
invoking possible earlier mis identification of severa l rapidly- 
evolving VLBI components (Marscherl ll996l Il999l) . another 
analysis based on a dense time coverage (15 observing epochs 
over 3 years) confirmed the existence of only subluminal appar- 
ent motion (< 0.3c), therefore implying a very sm all viewin; 
angle to the line of sight of 0.4° for the VLBI jet JPiner et 



ng 
al. 



1999). This analysis has recently been refined with an extended 
data span (28 epochs over 8 years), leading to a revised appar- 
ent speed of only 0. 1 + 0. 02c for the fastest VLBI jet component 
JPiner & Edwardsll2005h. While this val ue is smaller than the 
value found by Keller mann et alJ ( 120041) . who reported an av- 
erage component speed of 0.4c from 15 GHz monitoring over 
6 years, the two results do rule out the existence of superlumi- 
nal motion in the Mkn 421 jet. Interestingly, Mkn 421 is not a 
unique case and such low ap parent speeds have been found in 
the other TeV blazars as well (tPiner & Edwardsl2004l) . 

The apparent inconsistency between the low degree of 
relativistic beaming derived from the VLBI observations of 
TeV blazars and the high value predicted by the theory 
may be explained if the Doppler factor decreases along the 
jet as a result of either jet curvature or jet deceler ation 
JGeorganopoulos & Kazaiiasl 120031 iGhisellini et ail 120051) . It 
is also possible that the measured VLBI jet speed does not 
correspond to the speed of the actual underlying jet but in- 
stead to the speed of a perturbed pattern along the jet jZensua 
119971) . In this case, the radio core may still be efficiently 
boosted although there would be no evidence for relativistic 
beaming from apparent motions in the VLBI jet. Nevertheless, 
one would expect correlations between radio core properties 
and high-energy events if this hypothesis is correct. One in- 
dication of the existence of these correlations comes from 
the apparent connection between the epoch of emergence of 
new VLBI components and the occurrence of strong X-ray 
and y - ray flares in several y-ray emitt i ng AGN dWehrle et al] 
19931 IPohletalJ Il995t lUnwin et alJ Il997t btterbein et al 



1^0 ; irom ei ai .i iiy^ji junwin ei aij iiyy/i luti eroein ei ai. 
1998 iBritzen et all 1998: Bowerll998UKrichbaum et alll998l 



WehrlelH9 99: Marchenk o et al.ll2000l) . Furthermore, it appears 
that the y-ray flares detected by the EGRET detector onboard 
the Compton Gamma-Ray Observatory seem to occur during 
the rising phase of high-frequency radio outbursts, again sug- 
gesting a p ossible connection between radio and high-energy 
properties dValtaoia et al J 1 1 996l IValtaoia & Terasrantal fl996: 
Lahteenmaki et al. 200ojh This is especially true for Mkn 421, 
which occasionally shows multi-s pectral flares (from radi o 
to TeV energies), as reported bv iKatarzvnski et all (2003). 
Searching for additional clues along these lines with new multi- 
frequency data is important for investigating whether TeV 



blazars are indeed strongly-beamed sources, as presumed so 
far. 

This paper reports multi-frequency VLBI maps of Mkn 42 1 
obtained at three closely-spaced epochs (one month apart) 
in 1998. These observations were arranged during regular mon- 
itoring of the source at very-high y-ray energy (TeV level) by 
the CAT Cherenkov imaging telescope. In the following sec- 
tions, we present results of the VLBI and TeV observations, 
discuss the variability of the source, and investigate possible 
connections between VLBI properties and TeV activity. We 
also derive an upper limit for the projected size of the y-ray 
emission region based on the size of the compact radio core as 
measured from our high-resolution VLBI data. 

2. Observations and data reduction 

2.1. VLBI observations 

Multi-frequency (5.0 GHz, 8.4 GHz, 15.4 GHz, 22.2 GHz) 
polarization-sensitive VLBI observations of Mkn 421 were car- 
ried out at three epochs on March 4, March 28, and April 26, 
1998, using the Very Lo ng Baseline Array (VLBA) tele- 
scope (Napier et all 1 19941) of the National Radio Astronomy 
Observatory (NRAO) 1 . Dual-polarization (LCP and RCP) and 
two intermediate frequency channels (IFs), each 8 MHz wide, 
were recorded simultaneously at each band with scans of 3 min 
at 5 GHz and 8 GHz, 5 min at 15 GHz, and 6 min at 22 GHz. 
Mkn 421 was observed alternately with another TeV-emitting 
source (Mkn 501) at each of the four frequencies, along with 
interleaved calibrators for amplitude and polarization calibra- 
tion. The total observing time for Mkn 421 was roughly 0.7 hr 
at 5 GHz and 8 GHz, 1.2 hr at 15 GHz, and 1.4 hr at 22 GHz. 

The raw data were correlated with the VLBA correlator at 
the Array Operations Center in Socorro, New Mexico, and cali- 
brated with the NR AO Astronomica l Image Processing System 
in the standard way Juivestadl2000l) . The initial amplitude cal- 
ibration for each of the IFs and polarizations was accomplished 
using system temperature measurements taken during the ob- 
servations and NRAO-supplied gain curves. The instrumental 
polarizations for each antenna were determined from obser- 
vations of the unpolarized source OQ208 (1404+286), while 
the absolute polarization position angles were calibrated using 
observations of the compact polarized source 1823+568. For 
each epoch and frequency, the initial polarization position an- 
gles have been rotated so that those for 1823+568 match the 
integrated values reported in the University of Michigan Radio 
Astronomy Observatory data base 2 for our observing epochs. 
The rotation values derived with this scheme (equivalent to the 
R-L phase differences of the reference antenna) were found to 
agree within 2° for the data of March 4, 1998, and April 26, 
1998, at all four frequencies. This is expected if the R-L phase 
differences of the receivers are stable over this timescale; in 
fact, the stability of this phase difference has been demon- 
strated fOT_the_VI i BAantennas on timescales as long as two 
years (R eynolds et alJ2o"oih . It was intriguing, though, that the 

1 NRAO is operated by Associated Universities, Inc., under coop- 
erative agreement with the National Science Foundation. 

2 see http://www.astro.lsa.umich.edu/obs/radiotel/umrao.html 
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rotation values derived for the data of the intermediate epoch 
(March 28, 1998) were inconsistent with these two sets of val- 
ues. A closer look at our calibrator data for March 28, 1998, 
revealed afterwards that these data are indeed suspicious, prob- 
ably because the polarization properties of 1823+568 varied 
over the 12-hour span of our observations. Therefore, we cali- 
brated the absolute orientation of the polarized position angles 
for the data of March 28, 1998, using the mean of the rotation 
values determined for the other two epochs. 

The calibrated visibility data of Mkn 421 for each fre- 
quency band were processed using the Brandeis software to 
produce hybrid maps of the distribution of the total intensity / 
and linear polarization P 3 . Maps of the linear polarization were 
made by referencing the calibrated cross-hand fringes to the 
parallel-hand fringes using the antenna gains determined in 
the hybrid-mapping procedure, then Fourier-transforming the 
cross-hand fringes, and CLEANing. One by-product of this 
procedure is to register the / and P maps to within a small frac- 
tion of the beamwidth, so that corresponding / and P images 
may be directly superimposed. Subsequent to imaging, circular 
Gaussian models were fitted to the self-calibrated visibility data 
using a non-linear least-squares algorithm implemented in the 
Brandeis software. Such models are useful for estimating the 
spectrum of the individual VLBI components and for tracking 
temporal changes in the source structure, as discussed below. 

2.2. TeV measurements 

Very-high y-ray energy measurements on Mkn 421 were 
obtained using the ground-based Cherenkov telescope CAT 
("Cherenkov Array at Themis") located in the French Pyrenees. 
This telescope, which has a reflector area of 17.8 m 2 and a very 
high-definition camera with 546 pixels of 0.1° x 0.1°, detects 
the Cherenkov light emitted by secondary particles resulting 
from cascades of p rimary high-energy y rays and cosmic rays 
in the atmosphere jBarrauetalJ ri998). In this technique, air 
showers induced by y rays are efficiently separated from the 
background of c osmic-ray showers by the shape and direction 
of their images JLe Bohec et aljfl~998h . The energy threshold 
of about 250 GeV at zenith is achieved by using fast photo- 
tubes and fast coincidence electronics, with an integration time 
of 12 ns matching the short-lived Cherenkov flashes. 

Mkn 4 21 was the first ex tragalactic source detected at the 
TeV level (Pu nch et alJ ll992). so on after its iden tification as a 
y-ray source by EGRET in 1992 llLinetaljf l992). Monitoring 
by CAT was conducted from December 1996 (when this tele- 
scope began operating) until June 2003 (when it ceased oper- 
ations). During this period, Mkn 421 was observed each year 
from November to June (when the source is visible during night 
time at the CAT site) for about two weeks per month (centered 
on the new Moon) within technical and weather constraints. 
A detailed descript ion of the y-ray si gnal extraction and data 
analysis is given in Pir on et alJ (12001). Light curves have been 
produced for the integral flux above 250 GeV by assuming a 

3 P = pe 2 '* = mle 2 '*, where p = ml is the polarized intensity, m 
is the fractional linear polarization, and x ' s me position angle of the 
electric vector on the sky. 



differential y-ray spectrum index of -2.9 representative of all 
available spectral measurements for Mkn 421. 

3. Observational results 

3.1. VLBI morphology and variability 

In order to provide an overview of our results in compact 
form, we present our images of Mkn 42 1 at all four frequen- 
cies and all three epochs in Figs. [0 and |2] The image param- 
eters are listed in Table [2 The total-intensity morphology in 
our images is very similar to that observed in previously pub- 
lished maps at these freque ncies , such a s in the large series of 
image s presented by IPiner et all i 19991) and IPiner & Edwards! 
(2005). We identify the core and the three long-lived jet com - 
ponents C6, C5, a nd C4 first detecte d bv IPiner et all (fl999). 
The latest epoch of IPiner et alJ i 19991) was in December 1997, 
only three months before our observations, and the cross- 
identification between our components and theirs is straight- 
forward. Accordingly, we have retained this naming conven- 
tion for ease of comparison. The additional 5 GHz component 
label ed CO is a combi nation of components C3, C2, and CI 
from lPineretail Jl999l) . 

As noted above, in order to derive a description of the 
source structure that could be used to quantitatively compare 
our results at the different frequencies and epochs, we obtained 
circular Gaussian model fits to all the calibrated / and P visi- 
bility data. The resulting Gaussian model components are pre- 
sented in Table 13 The errors reported in this table are our best 
lcr error estimates, as derived by examining the x variations 
when forcing changes in the individual model components (er- 
rors of lcr correspond to an increase in the best-fit^- 2 by unity). 
In some cases these errors appear reasonable, but in some other 
cases they clearly underestimate the true errors. For example, 
reasonable errors in the separation from the core r are proba- 
bly no less than about l/20th of the beam size (about 0.08, 0.06, 
0.03, and 0.02 mas at 5, 8, 15, and 22 GHz, respectively), which 
is larger than the errors reported for C6 and C5 in Table [2] es- 
peci ally at the lowest freque ncies (5 and 8 GHz). 

IPiner & Edwards! (120051) conclude that the jet components 
C4, C5, and C6 were present over the entire time range spanned 
by their observations (from 1994 to 2002) and that they showed 
at most subluminal speeds of expansion from the VLBI core 
(linear least-squares-fit speeds of 0.089 +0.066c for C4, 0.095 + 
0.024c for C5, and 0.029 + 0.026c for C6). Our images are con- 
sistent with these results, and show no evidence of systematic 
motions in the locations of any of these jet components over 
the two-month time interval covered by our observations. Thus, 
our data exclude the possibility that the source structure actu- 
ally evolves much more rapidly than deduced so far, but that the 
time intervals between previous VLBI images were too small 
to track the jet-component evolution. 

Our measurements in Table[2]also do not show any conclu- 
sive evidence of variability in the fluxes or polarizations of the 
jet components over our three epochs, except for possible mild 
changes at the level of 10-20 mJy in the total flux of the jet 
component C6. In contrast, we have a strong indication of vari- 
ability in both the total flux and the polarized flux of the VLBI 
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Fig. 1. Total-intensity maps of Mkn 421 at four frequencies (from upper to lower panels: 5 GHz, 8.4 GHz, 15.4 GHz, and 
22.2 GHz) and three epochs (from left to right: March 4, March 28, and April 26, 1998). Image parameters are listed in Tabled 
The FWHM Gaussian restoring beam applied to the images is shown as a cross in the upper left-hand corner of each panel. 
Locations of the circular Gaussian components fitted to the visibility data at each frequency, as reported in Table|2] are marked 
in the left-hand panels. 
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Fig. 2. Polarized-intensity maps of Mkn 421 at four frequencies (from upper to lower panels: 5 GHz, 8.4 GHz, 15.4 GHz, and 
22.2 GHz) and three epochs (from left to right: March 4, March 28, and April 26, 1998). Image parameters are listed in Table^] 
The FWHM Gaussian restoring beam applied to the images is shown as a cross in the upper left-hand corner of each panel. 
Locations of the circular Gaussian components fitted to the visibility data at each frequency, as reported in Table|2] are marked 
in the left-hand panels. 
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Table 1. Parameters of the total-intensity and polarized-intensity VLBI images of Mkn 421 plotted in Figs.^&El 



Beam" Total-intensity maps Polarized-intensity maps 



Frequency 


Epoch 


a 


D 


4> 


reaK 


Lowest contour^ 


reaK 


Lowest contour^ 


(OHz) 




(mas) 


(mas) 


( ) 


(mjy/beam) 


("to ot peak) 


(mJy/beam) 


(70 ot peals) 


5.0 


i 


2.93 


1.92 


-3.3 


357 


0.18 


3.98 


35 




2 


2.93 


2.08 


-11.0 


324 


0.25 


5.95 


12 




3 


2.89 


1.77 


-6.4 


317 


0.25 


1.46 


55 


8.4 


1 


1.76 


1.15 


-3.3 


358 


0.18 


4.54 


16 




2 


1.79 


1.17 


-12.0 


308 


0.25 


4.29 


17 




3 


1.63 


1.07 


-6.9 


315 


0.25 


4.21 


22 


15.4 


1 


1.00 


0.67 


-2.3 


325 


0.35 


6.07 


28 




2 


0.95 


0.63 


-8.1 


296 


0.50 


7.92 


17 




3 


0.92 


0.58 


-4.7 


295 


0.35 


4.72 


34 


22.2 


1 


0.72 


0.46 


-0.8 


249 


0.50 


5.53 


32 




2 


0.82 


0.49 


-14.8 


255 


0.71 


4.18 


48 




3 


0.69 


0.40 


-5.5 


231 


0.71 


4.41 


48 



" The restoring beam is an elliptical Gaussian with FWHM major axis a and minor axis b, with the major axis in position angle if> (measured 
from North through East). 
b Successive contour levels are each a factor of the square root of two higher. 



core at all four frequencies over the roughly two months cov- 
ered by our three epochs. In particular, the total flux of the core 
decreased by 30-70 mJy (corresponding to a change of 10- 
20%) between March 4 and April 26 at all four frequencies, and 
this decrease was accompanied by a rotation of the core polar- 
ization angle x by 30-60° at all frequencies except 5 GHz (see 
Table|3J. Although the core^ was relatively constant at 5 GHz, 
the corresponding core polarized flux was probably variable. 
These variations are discussed in more detail below. 

3.2. TeV variability 

An extensive study of the CAT measurements acquired on 
Mkn 4 21 between 1996 and 2000 is found in IPiron et alJ 
J200ll) . As reported in this paper, the level of the TeV emission 
for Mkn 421 changed significantly from year to year. While at a 
low level in 1996-1997 (mean flux above 250 GeV: <D>250Gev = 
2.5 + 0.9 x 10 _11 cm _2 s _1 ), it was more than twice as high the 
next year, when our VLBI maps were obtained (O>250GeV = 
6.0 + 0.6 x lO^'cnrV 1 in 1997-1998). Later on, the very- 
high-energy emission decreased to a level comparable to that 
of 1996-1997 (a> >2 50Gev = 2.7 ± 0.9 x lO^WV 1 in 1998- 
1999), whereas intense activity started again in January 2000. 

The nightly-averaged integral flux above 250 GeV obtained 
by CAT between December 1997 and May 1998 is plotted 
in Fig. Also indicated in this figure are the epochs of our 
three VLBI observing runs. The most striking features in the 
light curve are the two y-ray flares detected in January and 
March 1998, the second of which occurred on March 25, just 
three days prior to our second VLBI run (Fig.|4}. Additionally, 
the overall TeV activity seems to show a regular increase from 
the first to third VLBI observing epochs, although such a trend 
is difficult to establish because of the sparseness of the TeV 
data, especially near the first VLBI run (due to poor weather 
conditions). Most noticeably, the level of the TeV emission 



near our third VLBI observing epoch (April 26) is significantly 
higher than the mean TeV flux value for the campaign (materi- 
alized by the horizontal dashed-line in Fig. [3}, while that near 
our second VLBI run (March 28) is roughly at the same level. 

4. Discussion 

4.1. Variability of the VLBI core 

As noted previously, both the total flux and the polarization 
characteristics of the VLBI core were variable in the time 
spanned by our observations. In this section, we first discuss 
the variability in the total flux and its frequency dependence 
and then examine the variability of the polarization flux and 
position angle based on the model-fitting results in Table|2] 

Figure|5]presents the spectra of the core and the jet compo- 
nents C6, C5, and C4. The core and C6 were detected at all four 
frequencies, while C5 and C4 were detected only at the lowest 
frequencies (5, 8, and 15 GHz for C5; 5 and 8 GHz for C4). The 
fluxes for C5 and for C4 at the three different epochs all agree 
within about 10 mJy or less and show no indication of variabil- 
ity. Those for C6 show differences up to 20 mJy, especially at 
5 and 8 GHz between our last two epochs. However, we cannot 
say for sure whether these are real or due to systematic devia- 
tions in the flux estimates. Overall, there is no convincing evi- 
dence of variability in the jet-component fluxes. The spectrum 
of C6 is consistent with this component being optically thin 
throughout the range of frequencies we have observed (spec- 
tral index of about -0.6), as is expected for jet emission. The 
spectra of C5 and C4 instead appear to be flatter (spectral index 
of about 0.0 for C5 and -0.1 for C4). 

In contrast, the core total flux is appreciably variable when 
compared to the constancy of the flux densities of the jet com- 
ponents C6, C5, and C4. As shown in Fig. [5] the core flux den- 
sity decreases at all four frequencies from our first to our third 
epochs. The total changes are 56, 68, 28, and 38 mJy at 5, 8, 15, 
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Table 2. Model-fitting results. 
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2.4 


4.5 


1.2 


-116.9 


6.7 
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0.05 
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" The components fitted to the complex visibility data are of circular Gaussian form with FWHM axis d, total intensity /, and polarized 
intensity p. They are separated from the (arbitrary) origin of the image by an amount r in position angle 6, which is the position angle 
(measured from North through East) of a line joining the components with the origin. The direction of the electric vector on the sky is given 
for each component by its position angle x (measured from North through East). 

* Combinations of polarizations for C6+C5. 
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Mkn 421 light curve in March 1998 
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Fig. 3. Nightly-averaged integral flux above 250 GeV for 
Mkn 421, as measured by CAT between December 1997 and 
May 1998. Arrows stand for 2<x upper-limits when no signal 
was detected. The horizontal dashed line shows the mean flux 
over the 5-month period. The epochs of our three VLBI runs 
(corresponding to MJD=50876, 50900, and 50929) are indi- 
cated by vertical thick dashed lines. 



and 22 GHz, which are clearly significant, even if some of the 
core flux errors are underestimated. The decrease in the core 
flux is seen at each frequency, arguing against the possibility 
that it is due to calibration errors, since the data for each fre- 
quency were calibrated independently. In addition, if the flux 
density decrease were associated with calibration errors, we 
might expect it to be more prominent at the higher frequen- 
cies, while the opposite is true. Figure |5] also shows that the 
spectrum of the core varies as the flux density decreases. The 
core spectrum at our first epoch (March 4, 1998) has a maxi- 
mum near 8 GHz and is presumably self-absorbed at lower fre- 
quencies. The spectrum at our second epoch three weeks later 
(March 28, 1998) is complex, but roughly speaking has become 
flatter. The core spectrum at our third epoch (April 26, 1998) 
is instead similar to the peaked spectrum observed at our first 
epoch, though the flux density at 15 GHz is somewhat higher 
than we would expect in this case. Thus, it seems that the core 
spectrum first became flatter, then started to change toward be- 
coming peaked again, with an overall flux density level appre- 
ciably below its initial level on March 4, 1998 (our first epoch). 

There is also some evidence from our data that the polar- 
ized flux of the core may have been variable at 5 and 8 GHz, 
though it is difficult to be sure of this given the correspond- 
ing errors. Figure [6] plots the polarized fluxes of the core and 
the jet components C6, C5, and C4 for the three successive 
epochs. As in the previous case, the jet components were de- 
tected only at the lowest frequencies (5, 8, and 15 GHz for C6; 
8 and 15 GHz for C5; 5 and 8 GHz for C4). The core polariza- 
tion was clearly detected at 5 GHz at our first epoch (estimated 
flux of 3.2 + 0.5 mJy), together with polarized flux associated 
with the jet components C6 and C4; however, the core polariza- 
tion was not detected above the noise level at our second epoch 
(estimated flux of 0.5 + 1 .0 mJy), though the polarizations from 
the jet components C6 and C4 were detected and were equal to 
their values at our first epoch to within the errors (see Fig.|6}. 



VLBI 



4.6± 0.9 



50890 50892 50894 50896 50898 50900 50902 

Modified Julian Day 

Fig. 4. Nightly-averaged integral flux above 250 GeV for 
Mkn 421, as measured by CAT during the 10-day period prior 
to our second VLBI observation on March 28, 1998 (indicated 
by a vertical thick dash line at MJD=50900). Arrows stand for 
2cr upper-limits when no signal was detected. The horizontal 
dashed line shows the mean flux over this 10-day period. 

The core polarization was only weakly detected at 5 GHz at 
our third epoch (estimated flux of 1.6 + 0.5 mJy). Overall, 
the fractional core polarization over the t hree epochs is 1% at 
most, which is half of that measured bv lPollacketai](l2003l) 
in a 5 GHz image from December 2000. The 8 GHz mea- 
surements suggest that the core polarized flux also decreased 
at this frequency between our first two epochs, but formally 
the decrease does not exceed the error estimates (difference 
of 1.6cr). At 15 and 22 GHz, the uncertainties are too large 
to draw any firm conclusion about the core variability, even 
though the core polarized flux appears to be somewhat weaker 
at the third epoch than at the first epoch. The fractional core 
polarization at these f requencies is about 2% and is consistent 
with that measured by Mars cher et alJ J2002I) in a 22 GHz im- 
age from August 1997, but it is smaller than the average value 
of 4.6% reported bv lPiner & Edwards! d2005h from 22 GHz data 
acquired at four epochs in 2001-2002. Overall, none of the jet 
components show signs of variability (see Fig.[6jl. 

Figure[7]plots the variability of the polarization position an- 
gle^ for the core and the jet components C6, C5, and C4, based 
on the data in Table |2 As shown in this figure, the core polar- 
ization position angle at 5 GHz is stable over our three epochs 
and the orientation of the electric vectors at this frequency is 
aligned with the direction of the inner VLBI jet (about -35°). 

This value for the polarization angle of the core differs by 70° 

I 1 1 1 

from that measured by Pollack et al. (2003) in December 2000, 

an indication of the existence of significant variability on the 
long term even though none is detected over the two-month 
span of our data. In contrast, the core polarization position an- 
gles at 8, 15, and 22 GHz changed substantially in the time be- 
tween our first and third epochs, whereas no such changes were 
found for the polarization angles of the jet components C6, C5, 
and C4. Again, this differs from the behavior observed in 2001- 
2002, where significant polarization angle variabilit y was de- 
tected for component C6 but not for the core dPiner & Edw ards 
2005). In all cases, the rotations for the core exceed the 3cr er- 
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Fig. 5. Spectra of the VLBI core and jet components C6, C5, and C4 for our three VLBI observing epochs (March 4, 1998, 
plotted in red, March 28, 1998, plotted in black, and April 26, 1998 plotted in green). The flux errors correspond to the formal 
lcr errors from the model-fitting results in Table|2] The same flux scale is used in each panel. 
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Fig. 6. Polarized flux of the VLBI core and jet components C6, C5, and C4 for our three VLBI observing epochs (March 4, 1998, 
plotted in red, March 28, 1998, plotted in black, and April 26, 1998 plotted in green). The flux errors correspond to the formal 
lcr errors from the model-fitting results in Tabled The same flux scale is used in each panel. 



rors in the polarization angles; the most dramatic rotation is 
through more than 60° at 15 GHz, corresponding to about 10cr 
(see Table The direction of the rotation is the same at all 
three frequencies, with the x values becoming more negative 
and showing a roughly linear dependence on time. Since the ro- 
tation is not larger at the lower frequencies, it cannot be associ- 
ated with variable Faraday rotation in the VLBI core. The^- ori- 
entation at 15 GHz at our first epoch (-10°) is nearly aligned 
with the direction of the inner VLBI jet (about -25°), and the 
X orientation at 22 GHz at our third epoch (-117°) is nearly 
perpendicular to this direction. However, the other core x ori- 
entations at 8, 15, and 22 GHz do not bear clear relationships to 
the direction of the inner jet. The observed core electric vectors 
at 15 and 22 GHz rotate nearly in synchrony, so that there is a 
large (^ 50°) and nearly constant offset between the x values 
at these two frequencies, whose origin is unclear; as discussed 



above, this offset does not appear to be associated with Faraday 
rotation in the VLBI core. 

The complex frequency dependence and variability of the 
core polarization position angles may be due to the presence of 
underlying components with different polarization properties 
within the region encompassed by the VLBI core. If such sub- 
components have different polarization position angles and/or 
vary on different timescales, a complex behavior arises natu- 
rally for the polarization angles of the VLBI core when the 
individual subcomponents are superimposed. Repeated mea- 
surements of the total integrated polarization for Mkn 421, 
which indicate two predominating polarization states over the 
time span of our observations (Fig. [HJ, give support to this sce- 
nario. As seen in Fig. [8] the two polarization states have elec- 
tric vectors that are either parallel to the VLBI jet (x = 150°) 
or perpendicular to the jet (x = 60°), hence corresponding 
to orthogonal magnetic fields if assuming optically-thin syn- 
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Fig. 7. Polarization angle variations in the VLBI core and jet components C6, C5, and C4 at 5 GHz (plotted in green), 8 GHz 
(plotted in red), 15 GHz (plotted in blue), and 22 GHz (plotted in black) over our three VLBI observing epochs (March 4, March 
28, and April 26, 1998). The polarization angle errors correspond to the formal lcr errors from the model-fitting results in Table[2] 
The same x scale is used in each panel. 



chrotron emission. In this scenario, the dominance of either 
state in the integrated polarization measurements would be the 
result of relative flux variations occurring on short timescales 
for subcomponents within the VLBI core. Alternately, the two 
polarization states may also result from variations in the opti- 
cal depth (thin or thick emission) of such subcomponents if the 
underlying magnetic field is solely perpendicular to the jet. 

4.2. VLBI properties and TeV activity 

Our VLBI data have revealed striking temporal changes over 
a few weeks for the compact radio core in Mkn 421. To our 
knowledge, it is the first time that such substantial variations on 
short timescales have been observed by VLBI in a BL Lac ob- 
ject (with the exception of intraday variability). As noted pre- 
viously, the complex temporal evolution for the VLBI core of 
Mkn 421 occurred at a time when the overall TeV y-ray activity 
for Mkn 421 was rising (Fig. O; furthermore, a TeV flare was 
detected just three days prior to our second VLBI run (Fig.|4}. 
Both results qualitatively argue in favor of a link between radio 
VLBI and TeV y-ray properties. 

Assuming that the y-ray radiation is nested in the compact 
VLBI core at 22 GHz (the highest available frequency), we can 
directly constrain the size of the high-energy emission region 
(projected on the plane of the sky) to be smaller than 0. 15 mas 
(our measured VLBI core FWHM at 22 GHz, see Table 01 or, 
equivalently, 0.1 pc (3 x 10 17 cm) if adopting a distance of 
125 Mpc for Mkn 421 (z = 0.031, H = 71 km s" 1 Mpc" 1 ). 
This upper limit for the projected size of the y-ray emit- 
ting zone is consistent with typical e stimates on the order of 
10 16 cm derived from variability data (Kataoka 2000) or from 
broadband spectrum modeling ( K atarzvnski et all2003l) . While 
our determination does not improve over such estimates, it is 
solely based on high-resolution radio-interferometric observa- 
tions and hence is independent of any theoretical model. 



In order to gain insight into the origin of the radiation for 
the compact radio core in Mkn 421, we compared its flux with 
predictions from high-energy emission models. Figure|5]shows 
the radio-to-TeV spectral energy distribution derived by fit- 
ting a one-zone synchrotron self-Compton (SSC) model to the 
data of a multi-frequency campaign on Mkn 42 1 that happene d 
to coincide with our third VLBI run ( Takahash i et alJ 2000). 
Superimposed on this theoretical curve, our measured VLBI 
core fluxes at 5, 8, 15, and 22 GHz (as reported in Tabled) 
are also plotted. It is striking that the flux of the VLBI core 
at the highest radio frequency (22 GHz) is perfectly consistent 
with the prediction of the model, whereas those at the lower 
radio frequencies (15, 8, and 5 GHz) depart from the SSC pre- 
dictions significantly. This finding is a strong indication that 
the compact core observed by VLBI at 22 GHz is dominated 
by the radio counterpart of the X-ray and y-ray SSC radiation. 
Conversely, the VLBI core at 15 GHz and below must encom- 
pass external radio-emitting regions located beyond the y-ray 
emitting zone, since the measured VLBI core fluxes at these 
frequencies are higher than the predictions of the SSC model. 

The above scenario, which interprets the 22 GHz radio 
emission from the VLBI core of Mkn 421 as the low-energy 
tail of a one-zone SSC radiation, naturally implies a connec- 
tion between radio variability and TeV activity. The existence 
of this connection is also supported by single-dish radio mon- 
itoring which revealed that some radio flares are correlated 
with TeV flares, with a simultaneous increa se in the radio and 
high- energy fluxes when such events occur jKatarzvnski et alJ 
2003J). In our case, however, we instead detect an overall de- 
crease in the radio flux of the VLBI core (Fig. [5} at a time when 
the very-high-energ v y-ray flux of the source tends to increase 
(Fig. 01. As noted by K atarzvnski et alJ J2003I) . self-absorption 
effects in the radio range can produce anticorrelation between 
flux variations observed at different energies within a SSC sce- 
nario. This may well be the case for our VLBI and TeV data, 
since the four observed radio frequencies are clearly within the 
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Fig. 8. Distribution of the position angle for the total inte- 
grated polarization at 4.8 GHz as observed by the University 
of Michigan within 100 days of our mid- VLBI epoch. Note 
two perpendicular states of polarization at about x — 150° and 
X = 60° in these measurements. 

self-absorbed part of the Mkn 421 spectrum (Fig. Under this 
assumption, the observed anticorrelated radio and very-high- 
energy variations would fit within the proposed SSC scenario. 

5. Conclusion 

Coordinated multi-frequency VLBI polarization and very- 
high-energy TeV observations of the y-ray blazar Mkn 421 
have provided evidence of variability in the compact VLBI core 
on a timescale of a few weeks and at a time when the overall 
TeV activity of the source was rising. Although we cannot be 
entirely sure about the connection between the two phenom- 
ena, our measurements are consistent with a scenario in which 
the VLBI core at 22 GHz represents the self-absorbed part of 
the low-energy emission induced by a one-zone synchrotron 
self-Compton model consistent with simultaneous high-energy 
data. At the lower radio frequencies (15, 8, and 5 GHz), 
the VLBI core encompasses additional emission unrelated to 
the SSC phenomenon and most probably originating from re- 
gions outside the y-ray emitting zone. Based on these observa- 
tions, we also derived an upper limit of 0.1 pc (3 x 10 17 cm) 
for the projected size of the y-ray emitting region, in agree- 
ment with previous estimates but based solely on our radio- 
interferometric data and hence free of any theoretical model. 

Further similar coordinated VLBI and TeV observations 
on Mkn 421 should be primarily targeted at improving the 
VLBI time coverage and angular resolution. Observations a 
few days apart over several weeks would be important to con- 
firm whether the VLBI core variability and TeV activity are 
indeed connected and whether there is any anticorrelation be- 
tween such variations. Increasing the VLBI angular resolution, 
e.g. by observing at 43 GHz with a transatlantic network or at 
86 GHz with the global millimeter VLBI array, would also be 
important to further constrain the location and size of the y-ray 
emission region in the framework of the proposed SSC sce- 
nario. By providing an angular resolution improved by a factor 
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Fig. 9. Multi-frequency spectral energy distribution from radio 
to TeV energies and VLBI core flux for Mkn 421. The con- 
tinuous line is a prediction based on the one-zone synch rotron 
self-Compton model derived by Takahash i et alJ (J2000) from 
the data of the multi-wavelength campaign conducted in late 
April 1998, while the dots represent our measured VLBI core 
fluxes at 5, 8, 15, and 22 GHz on April 26, 1998 (see Tabled). 
The inset shows an enlargement of the radio region. 



of up to 6, such high-frequency VLBI observations may per- 
haps even permit a direct measurement of the size of the SSC 
zone, which is on the order of 10 16 cm according to present 
models. VLBI monitoring of Mkn 421 should also be pursued 
on the long-term to refine jet component proper motions and 
possibly obtain more clues as to why there are no superluminal 
motions in this presumably strongly-beamed source. 
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